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ABSTRACT 

The blaclc-hole binary system LMC X-3 has been observed by virtually every 
X-ray mission since the inception of X-ray astronomy. Among the persistent 
sources, LMC X-3 is uniquely both habitually soft and highly variable. Using a 
fully relativistic accretion-disk model, we analyze hundreds of spectra collected 
during eight X-ray missions that span 26 years. For a selected sample of 391 
i?Ar£^ spectra we find that to within ^ 2 percent the inner radius of the accretion 
disk is constant over time and unaffected by source variability. Even considering 
an ensemble of eight X-ray missions, we find consistent values of the radius to 
within 4 — 6 percent. Our results provide strong evidence for the existence 
of a fixed inner-disk radius. The only reasonable inference is that this radius 
is closely associated with the general relativistic innermost stable circular orbit 
(ISCO). Our findings establish a firm foundation for the measurement of black 
hole spin. 

Subject headings: accretion, accretion disks — black hole physics — stars: indi- 
vidual (LMC X-3) — X-rays: binaries 



1. Introduction 



The X-ray binary LMC X-3 was discovered by Uhuru in 1971 (ILeong et al.lll97ll ). Ob- 
servations of its B3V optical counterpart revealed an orbital period of 1.7 days and a mass 
function of 2.3 ± 0.3 Mq. Because of its massive companion star, this established LMC 
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X 3 as a strong dynamical black-hole (BH) candidate (ICowley et al.l Il983l : iKuiper et al. 
19881 ). Subsequent X-ray observations spanning decades have r evealed a complex behavioral 
pattern that includes transitions b etween soft and har d states (jWilms et al.ll200l[ ) and long- 
term (> 100 d) variability cycles (ICowley et al.lll99l[ ). While by some metrics LMC X-3 is 
a nearly archetypal BH binary, its combined qualities of persistence and strong variability 
set it apart as unique. 

Among the black hole systems, LMC X-3 bridges the divide between low- mass X-ray 



binar ies powered by Roche-lobe overflow and wind-fed, high-mass X-ray binaries ( ISoria et al. 



20011 ). The former are transients, usually locked in a deep quiescent state, whereas the latter 
systems are persistently X-ray bright. Among the classical persistent BH sources (Cyg X-1, 
LMC X-1, and LMC X-3), LMC X-3 habitually shows the softest X-ray spectrum, reaches 
the highest luminosity, and exhibits the largest variations in intensity. 

Because of its persistence LMC X-3 has been observed by nearly every X-ray astron- 
omy mission. In this Letter, we apply our relativistic accretion disk model (kerrbb2; 



McClintock et al.l 120061 ) to essentially all available X-ray data in order to examine the pre- 
sumed constancy of the inner radius of the BH's accretion disk. We draw upon data collected 
by eight missions, with RXTE providing the lion's share. 

For thin accretion disks, recent MHD simulations provide support for id entifying the 



2008; 


Shafee et al. 


2008; 


Penna et al. 


2010 



a long history of theoretical and observational support (e.g., see Section 6 in lGou et al.ll2009l ). 



With this identi fication and the simple mon otonic relationship between -Risco and the BH 



spin parameter (IShapiro fc TeukolskyI 119831 ) . a measurement of Rm is equivale nt to a mea 



surement of the sp in of the BH. This i s the basis for both the continuum-fitting (jZhang et al. 



19971 ) and Fe-K (IFabian et al.l Il989l ) methods of measuring spin. In recent years, both 



methods have been us e d to estimate th e spins of stellar BHs (e.g.. IShafee et al. 



McClintock et al.l 12006 
2009h . 
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The mass of LMC X-3's BH primary is presently very uncertain (ICowleyl Il992l ). and 
we adopt a round value that is typical f or BH binaries of M = 10 M©. For the inclination 
we adopt the provisional value i = 67° ( iKuiper et al.lll988l ). Because of the uncertainties 
in both M and i, in this Letter we do not attempt to estimate the BH's spin. Rather, we 
assume a reasonable value for the mass and employ the X-ray continuum-fitting method 
in order to study the constancy of -Rin. We describe our data set comprised of hundreds 
of multi-mission spectra in Section |2l our analysis in Section |3l and present our results in 
Section HI In Section [5l we explore the systematics associated with our spectral model and 
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conclude by discussing our results in the context of ongoing studies of BH spin. 



Observations 



RXTE: The Rossi X-ray Timing Explorer {RXTE) is our workhorse instrument, provid- 
ing a total gross sample of 712 spectra. Individual spectra were defined by grouping all the 
archival pointed data from 1996 through 2009 into approximately half-day bins with 90% 
of exposure times ranging from 1-10 ks . We only use puls e-height spectra obtained by the 
best-calibrated PCA detector, PCU-2 ( iJahoda et al.ll2006l ). Count rates have been renor- 
malized to correct for detector dead time and a systemat ic error of 1% has been included 
to account for uncertainty in the instrumental response ( jJahoda et al.l 120061 ). These data 
have been analyzed from 2.55-25 keV over all reliable gain epochs (> epoch 2). Here and 
elsewhere, the analysis work has been performed using XSPEC v.l2.5.1o (recent enough that 
an early coding error in kerrbb has been fixecO; Arnaud 1996 ). 

EXOSAT: Seven observations from 1983-1984 were obtained via the HEASARC archiv^; 
only data from the ME instr ument are currently available. Spectra were extracted as de- 
scribed in [Treves et al.l ( 119881 ) and analyzed from 1-25 keV. The customary systematic error 
of 1% was included. 

Ginga: The LAC detector observed LMC X-3 on 18 occ asions during 1 9 87-19 90. To 
extract these spectra, we followed the procedures described in lEbisawa et al.l (119931 ). Each 
spectrum has been analyzed from 1.5-25 keV with a 1% systematic error included. 

ASCA: LMC X-3 was observed twice, once on UT 1993 September 22 and later on UT 
1995 April 14. We extracted and separately combined spectra from the two CIS and two 
SIS instruments. Data were calibrated relative to the GIS-2 detector and analyzed from 
0.8-9 keV(GIS) and 0.6-9 keV(SIS) using a 2% systematic uncertainty. 



BeppoSAX: Following the standard reduction guide ( iFiore et al.lll999l ). we have gener- 
ated spectra for the narrow-field instruments from each of the 23 available observations. We 
extracted spectra using 8' apertures in the imaging instruments and used a fixed rise-time 
threshold for the PDS. For each observation we employed all usable LECS, MFCS, and PDS 
data. Throughout, we adopted the standard inter-detector floating normalizations calibrated 
relative to the MFCS. Data were analyzed from 0.12-4 keV (LECS), 1.65-10 keV (MFCS) 
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and 15-80 keV (PDS). A 1% systematic error has been included. 



XMM-Newton: All photon-counting data were severely piled up and therefore rejected 
because of uncertainties in the flux calibration. We use the single available 19 ks timing- 
mode observation of LMC X-3 obtained on UTC 2000 November 25. Because of the large 
number of accumulated counts, ~ 2 x 10^, uncertainties in the response of the detector are 
dominant, and we therefore included a 3% systematic error and fitted over 0.5 — 10 keV. 
Reduction and processing has been performed using XMM SAS vQ.O.cjfl. 

Swift: The sole XRT windowed-timing mode observation of LMC X-3, taken on UTC 
2007 November 26, has been procured and analyzed following the procedures outlined in 
Capalbi et al.l (|2005[ ). We rejected all the photon-counting data because they suffer from 
extreme pile-up. Calibration version 11 files have been used for the data reduction. In 
consultation with the Swift Help Desk, we have included an extra model component to 
account for an instrumental artifact near the Si edge around 1.7 keV. Analysis has been 
conducted over 0.4-10 keV using a 1% systematic error. 

Suzaku: Two observations were made on Dec. 22 2008 and Dec. 21 2009 (UT). The 
Suzaku attitude calibration was improved using the AEattcor routing. We applied the 
appropriate reduction procedures for a bright point sourc^. Pile-up was kept well below ~ 
3% by excluding the innermost 10" and 30" for the 200 8 and 2009 ob s ervati ons, respectively. 
In all other respects, we have followed the methods of iKubota et al.l (|2010[ ). including using 
their energy intervals and adopting a 1% systemat ic error. A fi xed cross-normalization of 
1.16 is used between XIS and HXD-PIN detectors JMaedalboOsh . 



2.1. Flux Calibration 



Just as deducing the radius of a star from its spectrum requires knowledge of its lu- 
minosity, in order to estimate the inner radius of an accretion disk it is also necessary to 
determine its luminosity. However, the measurement of X-ray luminosity is problematic in 
X-ray astronomy because of the significant fiux-normalization differences, often ^10%, be- 
tween missions. We address this issue by using the power-law spectrum of the Crab Nebula 
as measured by lToor fc SewardI (119741 ): F = 2.1 and = 9.7 photons s~^ keV~^ at 1 keV. 



■^http: / /xnim.esac.esa.int / sas / 

''http: / /space. mit.edu/CXC / software /suzaku/ 

'"^ http://www.astro.isas. ac.jp/suzaku/analysis/xis/pileup/HowToCheckPileup_vl.pdf 
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For each mission considered herein (excepting Swift; see Table 1), we either rely on 
the Crab calibration performed by the instrument tea m, or we compu t e a co rrection to the 
effective area by comparing the spectrum predicted by iToor fc SewardI (Il974[ ) to parameters 
obtained by analyzing proximate, archival observations of the Crab. Toor & Seward normal- 
ization coefficients /ts and slope differences AFts are presented for each mission in Table 1. 
This table also summarizes for LMC X-3 the gross number of observations available from 
each mission, Nohs, as well as the number of observations that meet our selection criteria, 
iVsei (Section EH). 



3. Analysis 



At energies above ~ 5 — 10 keV, the spectra of BH binaries in all states show a contribu- 
tion from a power-law component. This power law is widely attributed to inverse-Compton 
scattering of thermal disk photons by hot coronal electrons. The power-law model we employ, 
SIMPL, gene rates this Compton component by upscattering seed photons from the thermal 
component (ISteiner et al.ll2009bl ). 



The thermal and principal component of our model is kerrbb2, a thin accretion disk 
model that includes all relativistic effects, self-irradia tion of the disk ("returning radiation" ), 
limb darkening, and the effects of spectral hardening (jhi et al. 2005 : McClintock et al. 2006 ). 
During analysis, this latter effect is handled on the fly via a look-up table of the spectral 
hardening factor / for a given value of the disk viscosity parameter a (we adopt a = 0.01 
as default). These tables were computed using bhspec, a second relativistic disk model 



flDavis et all 120061 : IPavis fc Hubenvl 120061 ). 



Our fit to the thermal component of the spectrum effectively determines the solid angle 
subtended by the accretion disk: Q = 7r(i?in/D)^cos i, where D is the distance and i is the 
inclination of the accretion disk with respe ct to the line of sig ht. For D we use the average 
distance to the LMC, D = 48.1 kpc (e.g., lOrosz et al.ll2009l ). while for inclination we use 



i = 67° (Section 1). Finally, we express i?in in dimensionless form, rm = Rin/{GM/c^) using 
M = 10 Mq (Section 1). We have recently shown that the choice of M, i and D, which 
effectively sets the absolute scale for Tin an d the luminos i ty, is q uite unimportant for testing 
the stability of (see Fig. 3 and text in ISteiner et al.ll2009al ). (These values are crucial, 
however, when it comes to estimating the spin of the black hole.) 

Using our adopted values of the source M, i, and D, our source model has four fit 
parameters: two for KERRBB2, Ri^ and the mass accretion rate M, and two for siMPL, 
the photon index F and fsc, which is the fraction of disk photons that get re-directed via 
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scattering into the power law. Our full model is tbabs(simpl®kerrbb2), where tbabs 
models the effects of photoelectric absorpti on; we fix its sole pa rameter: A^h = 4 x 10^°cm~^ 
(jPage et al.ll2003l ). using abundances from IWilms et al.l ( 120001 ) . For kerrbb2 we include 
limb darkening and returning radiation effects, set the torque at the inner boundary of the 
accretion disk to zero, and fix the normalization to unity. We use the efficient, up-scattering- 
only version of SIMPL, and in Section [5] we show that this choice is unimportant. 



3.1. Data Selection 



Our preliminary analysis of all the data showed that for many spectra the power- 
law index F was essentially unconstrained, even for the BeppoSAX, EXOSAT, Ginga, and 
RXTE missions, which have the requisite coverage to detect this component. This is be- 
cause the source is relatively faint (< 50 mCrab) and its Compton power-law component is 
generally very weak, showing a median normalization /sc ~ 0.3%. The extreme dominance 
of the thermal component in LMC X-3 makes it an ideal source for accretion-disk studies 
such as this. 

Restricting our census to the 134 RXTE spectra for which the photon index is measured 
to a precision better than ar = 0.5, we find a strong clustering of values in the range 
F fti 2 — 2.6. For our baseline model we fix F = 2.35 which matches the constant index 



derived from 22 deep RXTE pointings by lSmith et al.l (120071 ). and in Section [5] we show that 
our results depend very weakly on this choice for 2 < F < 3. 

Meanwhile, three missions, ASCA, Swift and XMM, have no sensitivity above E s^lO keV, 
and therefore only very loosely measure the power-law normalization parameter, /sc- At the 
same time, a self consistent and fruitful analysis of the thermal and Compton components 
requires that /sc be sensibly constrained. Therefore, and because the power law is generally 
so weak, we impose an additional data- select ion requirement, namely that for each fit /sc 
falls within the lower 95% span of the RXTE rank-ordered values. 

We further adopt a goodness-of-fit requirement, x^/v < 2, and a lower limit on the 
Eddington-scaled disk luminosity, /d = -^D/i^Edd > 0.05. This la tter criterion rem oves any 



hard state data in which the disk is likely truncated at r > r^n (e.g.. lEsin et al.lll997l ). Finally, 
i n consonance with the thin-disk model employed, we select only data for which Id < 0.3 
( iMcClintock et al.lbooeh . 



Table 1. Data and Instrument Summary 



Instrument N^^^ N^^i'^ /ts*" AFts'' Ref. 



RXTE (PCU-2) 


712 


391 (568) 


1.097 


0.010 




Suzaku (XISO) 


2 


2 (2) 


0.98 


-0.01 


1,2 


Swift (XRT) 


1 


1 (1) 


1.10 


-0.04 


3^^ 


XMM (MOS-1) 


1 


(1) 


1.00 


0.01 


4,5,6 


BeppoSAX (MECS) 


23 


2 (23) 


0.95 


0.00 


7 


ASCA (GIS-2) 


2 


2 (2) 


0.97 


-0.01 


8,9 


Ginga (LAC) 


18 


7(11) 


0.94 


-0.02 


10 


EXOSAT (ME) 


7 


6 (6) 


0.98 


0.00 


11 



References. — (1) ISerlemitsos et al.l l2007l : (2) 

htt p:/ /heasarc.g sfc.nasa.gov/docs/suzaku/pro p_to ols/suzaku_td/suzaku_td.html; 
(3) iGodet et al . 2009; (4) Guainazzi ct al. 20091; (5) [stuh lingcr ct al. 
2006 ; (6) private co mmunication w ith Ignacio de la C alle; ( 71 iFiore et al.l 



1999 
1989 



(8) Makishinmet_^ 



(iD lParmar fc Smith 



198 



19961; (9) lEbisawalll996l ; (10) [Turner et al 



'^Number of selected observations. Parentheses indicate the selection 
numbers when high luminosities (jj > 0.3 are allowed (see Fig.[2]l. 

^Jts is the ratio of the Crab normalization to that of Toor & Seward 
and AF-pg is the difference between photon indices. 

'^The Swift values are derived from a comparison between RXTE and 
Swift observations of 3C 273. 
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Fig. 1. — top: Accretion-disk luminosity in Eddington-scaled units (M = 10 Mq) versus 
time for all the data considered in this study (766 spectra). Red arrows show RXTE data 
which are off scale. Data in the unsh aded region satisfy our thin-disk selection criterion 
{H/R < 0.1, which implies Zd < 0.3; iMcClintock et al.l 120061 ). The dotted line indicates 
the lower luminosity threshold (5% LEdd) adopted in Section 13.11 bottom: Values of the 
dimensionless inner disk radius Tin are shown for thin-disk data in the top panel that meet all 
of our selection criteria (411 spectra; see Section [3]T]) . Despite large variations in luminosity, 
Tin remains constant to within ^ 4% over time. The median value for the RXTE data alone 
(rin = 3.77) is shown red dashed line. 
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4. Results 

The top panel of Figure [1] shows a 26-year record of the disk luminosity of LMC X-3, 
which is seen to vary by orders of magnitude. Two-thirds of the data meet our thin-disk 
selection criterion < 0.3. In the lower panel, we show the time history of the inner disk 
radius for just those data that meet all of our selection criteria (Section 13. ip . The radius 
is constant over the 26 years of monitoring to within ~ 2% for RXTE alone and ~ 4% 
considering all missions. 

Figure |2] explores the dependence of on luminosity. In this figure we include the 
high-luminosity data (/d > 0.3) that meet all of our other selection criteria (Section 13. ip . 
For Zd < 0.3 there is a gentle, nonlinear rise of Tin with luminosity. Especially visible in the 
RXTE data, this rise becomes prominent beyond Id ~ 0.25, above which there is a ~ 12% 
increase in rjn. No significant change in x^/v is associated with the apparent increase of 
Tin. We cannot say if this represents a real increase in at high luminosities or is simply 
an artifact of using the thin-disk model, which is expected to b e increasingly inaccurate at 
higher luminosities ( Penna et al. 2010 : Abramowicz et aPboiol ) at which a transition may 



occur to an advective slim-disk accretion mode. Interestingly, however, despite this rapid 
rise, we note that the RXTE data appear tightly clustered along a well-d efined curye. We 



approximate this dependence using a non-parametric curve-fit (LOWESS; IClevelandlll979l ) 
that allows us to detrend the data. We conclude that results from all eight missions, including 
the high-luminosity data, are in agreement with one another to within 6%. 



5. Discussion 

Figure [2] clearly demonstrates the limitations of the thin-disk model at high luminosities. 
We further illustrate this point in Figure [3] using LOWESS fits to the abundant RXTE data. 
We vary, one-at-a-time, the model components and parameters of our baseline model, group- 
ing these trials into four separate "families." In order of increasing importance, these families 
are (1) column density A^h, (2) power-law index F, (3) choice of power-law model, and (4) 
a. Figure [3] illustrates the changes introduced by adjusting each family of settings. 

We highlight two conclusions from Figured (1) Our results are relatively insensitive to 
all settings with the single exception of the choice of a- viscosity; the value a = 0.1 increases 
significantly the dependence of rin on luminosity. (2) The positive correlation between r\-a_ 
and luminosity is generally present for all families over the full range of luminosity, but it 
becomes prominent only above 0.2 — 0.3. 

Inspecting the families of curves in Figure 3 from top to bottom, one concludes the 
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Fig. 2. — The dimensionless inner-disk radius r-m versus luminosity for the filtered data 
(Section 13. ip and our baseline model. The vertical black line shows our adopted thin-disk 
upper limit, Id = 0.3. As in Figure [H the red-dashed line shows the RXTE average below 
this limit. 
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following: As the first two families show, our results are insensitive to the choice of A^h 
and only modestly sensitive to the choice of F. In modeling the Compton tail component 
(third family), one sees that our results are essentially identical whether one uses our base- 
l ine up-scattering-on ly model SIMPL = siMPL-1 or a a two-sided scattering model siMPL-2 



(jSteiner et al.ll2009al ). while the results obtained using the standard power law model POW- 



ERLAW differ only modestly (< 5%). 

The fourth family considers the primary setting for BHSPEC, the viscosity parameter a, 
used to compute spectral hardening (Section [3]). Here, we examine several distinct cases: our 
fiducial value, a = 0.01 (dotted), the value a = 0.1 (Section [3l dark blue), and alternative 
stress prescriptions q;md = 0.1 (orange) and ap = 0.1 (green). The parameter a typically 
refers to viscosity in the disk which is proportional to the total pressure at the disk midplane. 
However, other choices exist such as "beta disk" and "mean disk" models in which and 
aMD respectively describe viscosities whic h scale proportiona lly to the gas pressure or the 



geometric mean of gas and total pressures ( jPone fc Davisll2008l ). Both latter options produce 



spectral hardening values quite similar to those obtained for a = 0.01. In conclusion, only 
the second option, a = 0.1, has an important effect on our results. 

Our results indicate that the value of the inner disk radius rjn - and hence spin - is 
stable over decades, as is expected given the minute effects of accretion torques on a BH 
over such a time scale. We also confirm that rjn is nearly independent of luminosity provided 
that the disk is geometrically thin. The stability of r^n over time (for < 0.3) despite 
large fluctuations in the mass accretion rate provides strong evidence that and -Risco are 
closely associated, as we tacitly assume in measuring BH spin (Section [1]). 

The inter-mission consistency of our results (~ 4% below Id < 0.3 and 6% overall) is 
very important for future X-ray continuum measurements of BH spin: For some transient 
BH sources (e.g., A0620-00 and GRS 1009-45) only one or a few spectra are available in the 
data archives. Our results for LMC X-3 show that, as long as the power-law component is 
reliably measured, even a single, suitable spectrum can deliver an estimate of the disk inner 
radius accurate to several percent, and thereby a reliable measurement of spin. 



This research has made use of data obtained from the High Energy Astrophysics Science 
Archive Research Center (HEASARC), provided by NASA's Goddard Space Flight Center. 
We thank Shane Davis and Laura Brenneman for their insightful analyses of the manuscript 
and suggestions. JFS thanks Tim Oosterbroek for assisting with BeppoSAX LECS, and Ig- 
nacio de la Calle for advice on XMM. The authors thank Tomaso Belloni and Ken Ebisawa 
for contributing reduced archival data from EXOSAT, and Ginga and ASCA respectively. 
JFS was supported by the Smithsonian Institution Endowment Funds and JEM acknowl- 



- 12 - 



T — r 




Fig. 3. — Four families of models showing how our baseline results in Figure [2] are affected 
when a single model component or parameter is varied. The black dotted line drawn with 
each family of curves represents our fiducial model: A^h = 4 x lO^'^cm"^; F = 2.35; SIMPL-1; 
and a = 0.01. The horizontal dashed line for each family is set by the average value of r^n 
(see Figures [T] & 12]) , and each family is offset by 30% for clarity. Each curve represents a 
LOWESS curve fit to the RXTE data alone. Both axes are scaled logarithmically. 
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